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TECHNICAL NOTE NO. 1255 


THE EEESCT OF CC^EESSIBILITY ON THE GEONTH 
OF THE LAMINAE BOUNDARY LAYBE ON 
LOW-DEAG WINGS AI® BODIES 
By H. Julian Allen and Gerald E. NltzBerg 


SUMMARY 

The development of the laminar boundary layer in a compressi- 
ble fluid is considered. . Formulas are given for determinirog the 
boundary— layer thickness and the boundary— layer Eeynolds number, 
which la a measure of the boiuadary—layer stability, for airfoils 
and bodies of revolution. 

It is shoim that low drag coefficients can. be maintained to 
considerably larger Eeynolds numbers if these Reynolds numbers are 
associated with high ^fech numbers rather than low Mach numbers. 

The primary cause of the 3.noroase in boundary— layer stability with 
increasing Mach number is viscosity changes resulting from aero- 
dynamic heating. 


INTRODUCTION 

Experiments with a large number of low-drag airfoils have shown 
that as long as the transition from laminar to' turbillent flow at the 
sTirface occurs between the minimum pressure position and th^ trailing 
edge of the airfoil, the low— drag characteristics of these airfoils 
are maintained, but that as the transition point moves forward of 
the minimum presei.ure position the drag coefficient increases more or 
less markedly depending on the airfoil pressui’e distribution. 

It has been found that the boundary— layer Eeynolds number E 5 , 
based on the boundary— layer thickness and the local Velocity 
outside this layer, gives a fair measure of the stability of the 
boundary layer and, in conseq^uence, may be used as a criterion for 
determining the point at which transition to turbulent flow takes 
place. , , - 
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As noted in reference 1, the "best estimates of the critical 
value of Eg available at present were ‘obtained frcm fli^t tests 
of an NACA 35-215 airfoil sectidh ■which was tested as a. glove on 
the B— l8 airplane. Designating 



where V is the velocity outside the boundary layer; b, the 
distance from the surface of the airfoil to a point in the boundary 
layer where the velocity has reached 0.70? V; and v, the kinematic 
viscosity of the fluid, critical values of Eg between ^00 and 
9500 .were observed in these flight tests. 

Prom von Harman's momentum relation (reference 2, p. 107) it is 
evident that, in an incompressible fluid, if (l) the boundary layer 
on a given body is laminar from the stagnation point to any given 
point on the body and (2) the boundary— layer velocity profile is at 
all points along the surface of the same form when considered non- 
dimens ionally in terms of 5 and Y; then at the given point on 
the surface, the boTuidary— layer Heynolds number Eg is related to the 
body Eeynolds number E (baaed on body dimensions and the stream 
speed) in the form 

— - constant 
E 

If the constant is known for any body at the minimum pressure 
point it is possible to determine the body Eeynolds number, which is 
the upper limit of the range of the lovr-drag coefficients, for a 
given ■valtie of nearly incompressible flow this cons’tant 

may be e'valuated by the metliod of reference- 1. In those applications 
where the Mach number is not negligibly small, it is necessary to 
extend this- method to take account of the compressibility effects. 

Such an extension of this method is the subject of this paper. 


IHEOEY 

The growth of the laminar boundary layer in a compressible 
fluid may be conveniently studied by von TSS-rmixi'a momentum method. 
To this end, consider first the steady— s'fcate flow across the faces 
of an olemen'tal parallelepiped at the 8^q:faoe of a two-di.mensional 
body shown in figure 1. Let h, which is. chosen so as to be 
independent of' a, be the distance from the surface of ‘che body to 
a point in the boundary layer where the fluid shear has become 
negligibly small. 
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The several contributions to thq 3 ccaaponent of the change 
in momentuzn across the parallelepiped will now be considered in 
turn. 


The fluid entering the face normal to s per unit width 
introduces the moEientum 


I p dy 
^o 

while that removed at the opposite face is 
'■‘o 

and hence the change In thi,s contribution to the momentum is 

ds( / P u 2 dy ) - ( 1 ) 

No contribution occurs at the surface of the airfoil but at the 
parallel face an amount vT da is removed. Continuity requires 
that 


u2 dj + ds 


J:. 

da 




P U^ d3 


hence 


SXoyI = . . o Co.u). 
d y 3 b 


Py 


V = 


g 

ds 



p u dy 


and so this momentum contribution becomes 
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Since the flov is considered to be constant with time, the total 
change in the s component of momentum across the parallelepiped 
is given by the sum of eq.uations (l) and ( 2 ). 

The forces acting on the parallelepiped in the direction of s 
are the surface shear and the pressure difference between the 
surfaces normal to s. The shear force, using the established sign 
convention, is 


— T dS (3) 

and, if the boundary layer is thin, it has been shown (reference 2, 
p. 83) that the pressure variation with y is negligible, so that 
the pressure force Is 


-hfela (It) 

da , ' 

Wow V is sms.ll comcared to V. so that Bernoulli's equation for 
a compressible flow which Is cons'tant with respect to time may be 
written for the flow region outside the frictional influence of 
the boxmdary layer 

da da ■ 2 ds ’ 2 da 

It is convenient to rewrite this as 

^ (Pv^) , yg dpv ^ dV 

ds ds ds ds 


d (pyV^) , d (pyV) 
ds ds 


Moreover, since both py emd T are independent of y, then, for 
reasons which will be evident later, rewrite the last equaljicn as 


^ ^ _ 1 _d 
ds h ds 


p-V'^dy + .2. 


7 d 


h ds 'Jo 


I 


h 


Py7 dy 


'o 
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so tHat the expression for the pressure force hecames 


— h — ds 
ds 


ds 


ds 





ds 




<iy 


( 5 ) 


Finally, equating the change in the ■ s ccaaponent of momentiaa across 
the parallelepiped to the s directed forces on the parallelepiped 
the "momentum" relation for the tvo-dimensional flow of a compressi- 
ble fluid (i.e., with varying density) is 


T 



u®) dy - V 



u) 


<iy 


( 6 ) 


It has been observed in a number of experiments with conven- 
tional low-drag and high critical compressibility speed airfoils 
that the Blasiue— type profile is a good approximtion to the actual 
boimdary— layer profile over the forward region of the airfoil where 
the pressures are falling. An examination of the calculated 
boundary— layer profiles for a flat nlate at a number of Mach 
numbers (reference 3) Indicates that the form of the profile remains 
close to the Blaaius type for subsonic flows. As a consequence it 
seems reasonable to assume, as is done in the analysis to follow, 
that the boundarj’- layer over the surface of conventional low-drag 
and hi^ critical speed airfoils will remain of the Blasius type 
throu^out the subsonic speed range. 


For adiabatic conditions, the local temperature and density 
outside the boundary layer are, respectively. 


^v 






Pv " Po 





J 


and to the order of 


M 


the local density is 


Pv ~ Po 



( 7 ) 
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where the Buhscript o denotes conditions in the free stream and 
the subscript t denotes conditions Just outside the boundary layer 
at any point s along the airfoil, where the velocity is V; M is 
the free-stream Mach number, and y = Cp/cy, the ratio of specific 
heats , 

Since the- pressvire is transmitted unchan^d throu^ the boundary 
layer, it follows from the law of Boyle and Charles that the density 
at any point y within the boundary layer is related to the local 
temperature by 


Further, it is shown in reference 3- that for a flat plate the 
temperature variation within the boundary layer, for the Prandtl 
number equal to unity, is given by 


T 


= T. 


+ 



T. 


\ 


1 

. ^oj 


2 n 


The Prandtl number is denoted by 


Pr 


CpM 

k 


wherein, for the fluid, 

H the absolute viscosity coefficient 
Cp the specific heat at constant pressure 
and 

k the thermal conductivity 

For air, the Prandtl number is less than unity (at standard condi- 
tion Pr for air is O.733) but it is not expected tMt the form of 
the temperature •variation as given in reference 3 will, for air, bo 
seriously in error. For- the airfoil it seems correct then to assume 
the temperature vepriation to be of the same form 
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T 


T = T, 




— T ^ 
/ 


1 - 




( 9 ) 


Moreover, the results of tests with a circiaar cylinder (reference li-) 
have indicated that the surface temperature may he given with 
reasonable accuracy by 


T = 
■^u=o 


1 




( 10 ) 


Finally, from the relations of equv.tions (t), ( 8 ), ( 9 ), and (lO), to 
the order 1 ^, it may be found that 


P = Po ^ 




L 


1 + 


’ey 




The surface unit shear is gj.ven by 


T = 


du'^, 


dyy 

•T=o 


Experiment has shown that u varies as the absolute temperatur| 
to the 0.76 power and from aquations (t) and (lO) to the order of 



(7-1) 


arc--* 



and so, to the present order of approximation 


( 12 ) 


T = 1^0 





0. 76 (7-1) if r/x 

2 L 



(13) 


Using the density relations of equations (7) and (ll)^ the 
value of T given by equation (I3), and assuming the Blasius 
variation of u /7 with y/h in Cie momentum equation ( 6 ), it was 
found that, to the order of ly?, the boundary-layer thickness & 
is given by 
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where _ 

c chord of the airfoil 

Vq the kinematic viscosity in the free stream 

Yx the velocity outside the boundary layer at the chordwise stationj 
Si/c for which the boimdary layer is being computed 

6 the boundary -layer thickness^ which is considered in this analysis 
to be the distance from the surface of the airfoil to a point 
in the boundary layer where the ratio of the_local velocity 
to the velocity outside th:- boundary layer is O.7O7 

To employ E5 as a criterion -for determining the stability of 
the laminar boundary layer, accomt must be taken of the fact that 
because of aarodynamic heating, the kinematic viscosity varies 
throu^out the boundary layer. The value of V used in calculating 
the boundary— layer Eeynolds number shoulcUbe that characteristic of 
the point in the boundary layer at which instability initiates. The 
theoretical analysis does not indicate the location of this point, 
but experiment indicates (as will be discussed later), that insta- 
bility initiates near the inside of the boundary layer. The kine- 
matic viscosity is given by 



but equation (8) applies throu^ the boun^ry layer so 
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and using tha relations of equations (T) arid (dO) 


‘\i=o \ . 
Vo'/ 


1 + ^ M® Pr^y^T 

. «= \.''o/ 


n3L.76 


r / \ 211.76 

|_i J 


A - Z-ri 

I 2 


[©■-ir- •;> 


( 17 ) 


If 


_ ^UgQ _ \ / _JJ_ \ / \ 

Rc " <^^"0 V'^A^O / V’’u-o/ 


( 18 ) 


then with the vaJLue of 8 given by equation (ill-) 
r or/V\2 113.48 

|l-0.2Q^|(-) .-iJI ^ 
fl40.17ia^X /'li Y*^'" - ^ 

L Ao/'J Aoy 




1+1<? 


0.6T(;iV_o.35]} 


^Xj-^'"i(s/c) -o.Wif d(s/o) I ( 19 ) 


/tj-Yo.iv 


s/c) j 


For moderate Mach, numbers equation ( 19 ) may be approximated by 




0.41}- ^ ' f X d(s/c) 

V r 


( 20 ) 
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For the ccajrpreasihle fluid boundary 
tion, the momentum relation may he found 


layer of a body of revolu— 
to be 




r 


dy - V 



and under the prerious assumptions as tc temperature, and density 
variation and the shape of the boundary— layer velocity profile, it 
may be shown tliat to the order of 


&2 


Vi 


Z' V \ \ ® r f" 

(S fe) i”{ 


1 + M£ 


r 


r 


0.67(y-j)‘ 


- 0.35 


} 


° ( 22 ) 


■and 

% . \^v Yiy 


5 . 3 0 . 63 -0 • 3 ^^ 


Si/L 


° ( 23 ) 


where . 

ri radius of the body at si 

r radius of the body at s where the velocity is T 

L length of the body 



IIACA TN No. 1255 


11 


% = ^oL/Vo 

and the remaining sysibola arc as ^.-reviously designated. 

To apply the eq^uatione, the velocity distribution at the Mach 
number M must be ascertained. When the experimental pressure 
coefficient P distribution is known at the desired 1-iach number, 
the distribution of 'V,/7 q may be found using Bernoulli’s equation 
for a ccmpressible fluid. Por air this equation is ' 

f _l''f = 1 + 1 - [I-hO.7025 

^Vo/ C.2025 

Values obtained frcaa this equation are given in table I. 

For two-dimensional flow, whex-e the pressure coefficient 
distribution is knovm for M = 0, ihat for the desired Mach numbex- 
may be calcudatod using the von Karman— Tsien equation (reference 5) 


P « 


?M=0 






( 25 ) 


2 (1 +Vi-J-!F) 

Values obtained from this equation are givexx in table II . 


mSOTSSIOH fil’ih CONCLUSIONS 

An investigation of the boundary— layer thickness at a point 
55 percent of the chord behind the leading edge on the upper surface 
of an RACA 66, 2— il-20 airfoil at several Mach numbers was .conducted 
in the l6-foot wind tunnel at the Ames Aeronautical Laboratory. 

Using the measui-ed pressure distributions at the same Mach numbers, 
the boundary— layer thickness was calculated by equation ( 111 -) which 
considers effects of compx'’essibility and aerodynamic heating, and 
by the corresponding equation of reference 1 which neglects these 
effects. The calculated variation of boundary— layer thickness with 
bfe-ch number as determined from these equations and the several 
experimentally measured values are shovm on figixre 2, That the 
theoretical variation of 6 is valid is indicated by the close 
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agreement between the calculated ■Values obtained frcm eq,viatlon (ih) 
end the eqtperimen-tal results. 

Aa noted previously, in order that Bg may b© used as a cri- 
terion for -the stability of the bovmdory layer, it is essential to 
determine where, within, the bovindary layer, transition to turbulent 
flow initiates. In experimeii'tal investigations wi-th flat plates 
(reference 6), it was fotuid that slow fluctuations of flow occur 
within the boundary layer "thou^ they are not apparent near the 
outside of the layer. Jones (reference 7) ob-talned experlmen’fcal 
results substantiating "these data and suggested that the phencoena 
of transition to turbulent flow may be the direct result of inter- 
mittent ins-tabillty due to transient separation of the flow from 
the su2face. If this is true> then transition must Initiate near 
the Inside of the boundary layer. 

The experimental data a-milable showing the effects of ccstt- 
pressibility and aerodynamic heating Indicate that -toansition does 
arise near the surface. These da-fca were ob-fcained With four WACA 27-S12 
airfoil models <fef differeht chords by measuring -the maximum Eeynolda 
number for which low dreig was maln-fcained. It was found -that the 
values of this critical Eeynolds number .for , the smaller chord airfoils, 
which req.ulred higher ^fc,ch numbers "than •bhe larger chord airfoils to 
reach a given Eeynolds number,' were much higher than those for ■fche 
larger chord airfoils. The variation of critical Eeynolds number with 
Mach ntember computed from equation (I9) and experimental measurements 
for the NACA 27-^12 airfoils are shown in figure 3» It is seen that 
when the boundary^-layer Eeynolds number is based on the inside vis- 
cosity and the low Mach number experimen'fcal points ore fitted to the 
theoretical curve, the calculated effect of ocaapresslbility is in 
agreement with experiment. In -this connection, a calculation was 
made to determine the ratio ^®i^g "the viscosity coire— 

spending to conditions outside the ooundary. The calculation pre- 
dicted a decreasing •value of wl'th Mach nu^er, which is 

con-trary to "the observed fact. 

Eecently an experimen'tol investigation was mode of the effect 
of heating the surface of a low-drag airfoil by heating elements 
placed vi-thln -the wing. It was found, when heat was so applied as 
to main'tain ■fche entire surface of the airfoil over which the laminar 
flow occurred at a constant temperature increment above the tempera— 
txire of the ambient stream, -that the boundary layer yfaa destabilized 
so that the critical Eeynolds number "was decreased, Wi"fch "the same 
temperature increment but with only "the sxjrface in "the immedlato 
vicinity of "the transition region heated, the des"fcabllizlng effect 
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on the "boundary layer was even more marked. These results, in 
contrast to the results of the eocperiments with the NACA 27—212 
airfoils previously aJJLuded to, would indicate ■ that viscosity con- 
siderations alone are not stifficient to explain the effect of heat 
on the stability of the boundary layer unless, as in the cases of 
aerodynamic heating and no heating, the temperature gradient in the 
boundary layer at the surface is zero. ■ In the heated airfoil inves- 
tigation, heating only the leading-edge section was expected to 
produce a temperature variation in the boundary layer at transition 
similar to that obtained in aerodynamic heating. This was attempted 
but the heat so transfe3rred was insxifficient to materially influence 
the temperature variation -over and above that occurring naturally 
at the Ifech numbers of the tests so that no further conclusions could 
bo drawn. • ' 

An exact solution of the boundary-layer manentum equation can 
be obtained for compressible fl6w over a flat -plate. In figure h, 
the exact theoretical variation of the flat— plate critical Reynolds 
number is compared with values obtairad from equations (I 9 ) and (20) . 
It is seen from this figure that, even at supersonic Mach numbers, 
equation (I 9 ) is in good agreement with the exact solution. Con- 
siderable caution should bo used in applying the analysis of the 
present report at such large Mach numbers, ■ The boundary— layer 
thickness computed by equation (l4) should bo reasonably accurate 
oven at somewhat supersonic Mach numbers j but there is only scant 
experimental basis for the assumption that the same value of 
determines the limit of stability of the laminar boundary layer at 
sonic velocities and at low speeds. In consequence, it is con- 
sidered that the equations for stability developed in this report 
should only be used over the range of Mach numbers for which 
equations (I 9 ) and (20) are essentially in agreement. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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TABIS I 





Yeluos 

offX 

V^c 

~/rr 

V 



1"^. M 
1 

0.1 

0.2 

0*3 

0.4 

' 

0.5 

t 

0.6 

0.7 

1 J 

0.8 

1 0.7 

o.oon 

0.0043 

0.0097 '0.0168 

0.0255 

0.0358 

0. 0454 i 6. 0579 

' .6 

.0007 

.0028 

.0063 

.0109 

.0164 

.6232 

.0298! 

.03S1 

i *5 

.0004 

.0017 

.0039 

.C069 

.0105 

.0145 

.0194 

.0248 j 

1 

.0002 

.0010 

.0023 

.0041 

.0062 

.0087 

.0U7i 

.0149 

! .3 

.0001 

.0006 

.0012 

.0020 

.0(^1 

.0047 

■.00611 

.0077 I 

' .2 

.0000 

.0003 

.0006 

.0009 

.0013 

.0020 

.0027 

.0033 

1 .1 ^0 

.0001 

.0002 

.0004 

.0006 

.0008 

.0010' 

.0012 

: 0 

0 

0 

0 

0 

0 

0 

0 

0 

i -0.1 

0 

0.0000 

0.0002 

O.OOCQ 

0.0004 

0.0006 

0.0008 

0.0010 

; -.2 

0.0000 

.0001 

.0004 

■.0009 

.0011 

.0016 

■.0023 1 

.0031 

! -.3 

.0001 

.0003 

.0009 

.0016 

,0025 

.0037 

. 0052 ■ 

.0069 j 

i -A 

.0002 

.0006 

.0016 

.0027 

.0044, 

.0065 

.0089 

.0121 1 

: -.5 

.0003 

.0009 

.0023 

.0042 

.0069' 

.0100 

.0137, 


i -.6 

.0004 

. 0013 

.0032. 

.0059 

.00961 

.0141; 

.0195^ 

■ 1 


.0005 

.0018 

.0043' 

.0078 

.0127; 

.0187 

— -p -1.‘ 

. — ! 

! -.8 

.0007 

.0023 

.0055: 

.0100 

.0161; 

.0242} 

» .i. ■ 

_ 

! -.9 

.0008 

,0029 

.0068: 

.0124 

.0200- 

.0304> 



i -1.0 

.0009 

.0035 

.0083; 

.0151 

,0245! 

1 

— — — j 

— 

— 

1 -1.1 

.0011 

.0042 

. 0099' 

.0179 

.0292: 


. 

— . . 

, -1.2 

.0013 

.0050 

.0115j 

.0210 

.0340' 

1 

.! 

— 

1 
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Val^•0B of (P — Pm=o) 



0.1 ■ 

T 

1 0.2 
i 

0.3 

OA 

-O.5 

0.6 

B 

0.8 

0.7 

0.0022 

*0.0094 

0.0220 

0.0428 

0, 0660 

0.IG32 

0.1550 

0.234 

•6 

.0021 

j .0086 

.0201 

. .6392 

.06IQ 

.09-741 .1475 

.226 

•5 

i0019 

1 .0077 

■-.0179 

.0345 

.0556 


.214 

.h 

.0016 

1 .0064 

.0150 

.0297 

.0479 

.0758 

.1203 

.184 

•3 

.0012 

. 0052 

.0119 

.0231 

.0383 

.0598 

.09S9 

• .150 1 

.2 

.0008i ,.0037 


.0162 

.0273 

.0420 

. .0727 

.110 ; 

.1 

.0004 

.0020 


.00S2 

.0142 

.0221 


.0613 ’ 

0 

0 

0 

0 

0 

0 

0 

0 

0 : 

-0.1 

-0.0005 

-0.0022 

-0.0049 

-0.0098 

-0.0160 

-0.0271 

-0. 0450 

-0.0724 I 

-.2 

-.0011 

-.0046 

-,0104 

-.0201 

-.0346 

-.0562 

-.0970 

-157 i 

-.3 

-.0017 

■ -.0078 

-.016^1- 

-.0312 

-.0542 

-.0895 

-.1540 

-.2555 i 

-.4 

-.0024 

-.0102 

-.0230 

0438 

-■.0772 

-.1264 

-.2145 


-.5 

-.0032 

-,0130 

-.0305 

-.0571 


-.166? 

-.2776 

^ » 

—.6 

-.0040 

-.0161 

-.0378 

-.0706 

-.1273 

-.2080 

— 3440 

^ ^ » 

-•7 

-.0049 

-.0194 

-.0461 

-.0859 

-.1559 

-.2570 


^ ^ 

-.8 

-.0053 

-.0229 

-.0546 

-.1021 

-.1868 

-.30^ 


^ ^ 

-.9 

-. 0067 

-.0269 

-.0642 

-.1199 


-.3645 

^ ^ ^ 

M* VSk 

-1.0 

-.00761 

-. 0312 

-. 0742 

-.1382 

-.2515 

— 


^ 4^ 

-1.1 

-. 0086 j 

-.03491 

— 0855, 

-.1579 

-.2865 


» — — 

— 

-1.2 

-.0098 

-.03961 

-.09711 

-.1783 

-.3245 

" -J 


— — 

— 
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Fig. 2 NACA TN No. 1255 
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Maob number, M 

Figure 3.- Tbe effect of oompresBlbility on the critical 
Reynolds number ratio of the NACA 87-S12 

airfoil. 
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Fig. 4 
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Figure 4.- The effect of compressibility on the critical 
Reynolds number ratio of a flat plate. 



